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We study the stability and maximum mass of hybrid stars, assuming a generic quark matter equation 
of state with a single first-order phase transition between nuclear and quark matter, and a sharp 
interface between the quark matter core and nuclear matter mantle in a neutron star. For standard 
nuclear matter equations of state, we find that the mass-radius relation contains a stable hybrid branch, 
connected to the nuclear matter star branch, if the energy density discontinuity at the nuclear-quark 
transition is less than a critical value, which is typically between 60% and 80% of the nuclear matter 
energy density at the transition. 

Extending the quark matter equation of state to higher densities by assuming it has a density- 
independent speed of sound, we find that, as has been noticed before, there can be a disconnected 
branch of hybrid stars (sometimes called "third family" of stars). For typical nuclear matter equations 
of state, this branch exists if the nuclear matter density at the transition is less than a critical value 
which is about two to four times nuclear saturation density. We calculate the maximum hybrid star 
mass as a function of the nuclear matter density and the parameters of the quark matter EoS, and find 
that, for apparently reasonable values of the quark matter parameters, hybrid stars with mass above 
2 Mq can exist. 

PACS numbers: 25.75.Nq, 26.60.-c, 97.60.Jd, 



I. INTRODUCTION 

There has been much work on a plausible form for the 
equation of state (EoS) of nuclear matter at densities 
above nuclear density, using models of the nuclear force 
that are constrained by existing scattering data (see, for 
example, [IH3])- However, we remain almost completely 
ignorant of the quark matter EoS at the low tempera- 
tures and the range of densities that are of relevance to 
neutron stars. This is because cold quark matter cannot 
be created in laboratories, and numerical studies using 
the known strong interaction Lagrangian are stymied by 
the sign problem (see, for example, [1HZ])- 

In this work we therefore assume a generic quark mat- 
ter equation of state, and see in what ways it might be 
constrained by measurements of the mass and radii of 
compact stars. We assume that there is a first-order 
phase transition between nuclear and quark matter, and 
that the surface tension of the interface is high enough 
to ensure that the transition occurs at a sharp interface 
(Maxwell construction) not via a mixed phase (Gibbs 
construction) . This is a possible scenario, given the un- 
certainties in the value of the surface tension [8 , 9J . (For 
analysis of generic equations of state that continuously 
interpolate between the phases to model mixing or per- 
colation, see Refs. |10l ITT] .) 

We then answer the following questions: 

1. What is a sufficient condition for a completely 
generic quark matter equation of state to yield a 
stable branch of hybrid stars? 

2. If the quark matter EoS is not completely generic 



but the only assumption we make about it is that 
it has a density-independent speed of sound (like a 
classical ideal gas [29]), what form does the hybrid 
star branch take, and how heavy can those stars 
be? 

To address these questions we need to assume a nu- 
clear matter EoS, but for the first question we only need 
to know some generic features of the quark matter EoS, 
namely the pressure p — p crlt at which the transition from 
nuclear matter to quark matter occurs, and the disconti- 
nuity in energy density Ae at the transition (Fig.[TJ). Our 
answer to the first question will therefore depend on p cr it 
and Ae but it will be model-independent, requiring no 
additional assumptions about the quark matter equation 
of state. 

To address the second question, we will need the quark 
matter equation of state over a range of pressures above 
p cr i t . In this paper we assume, that the quark matter has 
a classical-ideal-gas equation of state, for which the speed 
of sound Cq M remains constant as the pressure varies 
from p cr it up to the central pressure of the maximum 
mass star. The equation of state for dense matter is 
therefore approximated as (see Fig. [IJ 
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where £nm(p) is the nuclear matter equation of state. In 
Appendix [X] we describe the thermodynamically consis- 
tent parameterization of the EoS that we used for quark 
matter. A similar generic parameterization is used in 
Ref. [15], which also considered the possibility of two 
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EoS 


max mass 


radius at M = 1.4 M Q 


L 


NL3 


2.77 M Q 


14.92 km 


118 MeV 


HLPS 


2.15 M 


10.88 km 


33 McV 



TABLE I: Properties of the NL3 and HLPS equations of state. 
L characterizes the density-dependence of the symmetry energy 
(see text). NL3 is an example of a stiff EoS, HLPS is an example 
of a softer one at density n < 4no . 



first-order transitions involving two different phases of 
quark matter. 
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FIG. 1: Equation of state e(p) for dense matter. The quark 
matter EoS is specified by the critical pressure p CIlt , the energy 
density discontinuity Ae, and the speed of sound in quark mat- 
ter Cq M (assumed density-independent). Below p cr it, the NL3 
EoS (curved line) is shown for illustration. 



no, there are many proposals in the literature. For il- 
lustrative purposes, we use two examples: a relativistic 
mean field model labelled NL3 [TB] and a non-relativistic 
potential model labelled HLPS, corresponding to "EoSl" 
in Ref. [SJ. Some of their properties are summarized in 
Table [T] where L is related to the derivative of the sym- 
metry energy S2 with respect to density at the nuclear 
saturation density, L = 3no(dS2/dri)\ no . We see that 
NL3 is an example of a stiff EoS, with a high value of 
L, producing stars with a high maximum mass and large 
radius. HLPS is a softer equation of state at n < 4tiq. 
There is some evidence favoring a soft EoS for nuclear 
matter: in Ref. [T7], values in the range L = 40 to 60 MeV 
were favored from an analysis of constraints imposed by 
available laboratory and neutron star data. Using data 
from X-ray bursts, Ref. [18 finds the surface to volume 
symmetry energy ratio S s /S v ~ 1.5 ±0.3 (See after their 
Eq. (43) and Table 4), which corresponds to L in the 
range 22 ± 4 MeV (using Eq. (7) in Ref. (17]). 

• Nuclear/ Quark transition. The nuclear matter to quark 
matter transition occurs at pressure p cr n . We will usually 
specify its position in terms of the baryon density n cr i t 
of nuclear matter at the transition, or the energy density 
Edit of nuclear matter at the transition. For a given nu- 
clear matter EoS, one can use any of these quantities to 
specify the location of the transition. We will typically 
give 7i cr it in units of no = 0.16 fm~ 3 . 

• Organization of paper. In Sec. ^ we discuss the criteria 
for stable hybrid stars to exist, as a function of the nu- 
clear matter equation of state and the parameters of our 
generic quark matter equation of state. We present re- 
sults for a hard and a soft nuclear matter EoS. In Sec. IIHI 
we discuss the maximum mass that such hybrid stars can 
achieve. Sec. |IV| gives a summary and conclusions. 



• Quark Matter EoS. The assumption that quark mat- 
ter has a density-independent speed of sound is rea- 
sonably consistent with some well-known quark matter 
equations of state. For some Nambu-Jona-Lasinio mod- 
els, the EoS fits Eq. (JT|) almost exactly [IS]. In addi- 
tion, the perturbative quark matter EoS [13 has roughly 
density-independent Cq M , with a value around 0.2 to 0.3, 
above the transition from nuclear matter (see Fig. 9 of 
Ref. |14|L In the quartic polynomial parameterization 
|15| . varying the coefficient 0,2 between ±(150MeV) 2 , and 
the coefficient 04 between 0.6 and 1, and keeping n CI i t /no 
above 1.5 (no = 0.16fm -3 is the nuclear saturation den- 
sity), one finds that Cq M is always between 0.3 and 0.36. 

In this paper we study the maximum hybrid star mass 
for a range of values of Cq M , from 1/3 (characteristic of 
very weakly interacting massless quarks) to 1 (the maxi- 
mum value consistent with causality) . We hope that this 
will give us a reasonable idea of the likely range of out- 
comes for realistic quark matter in which Cq M could vary 
with density. 

• Nuclear Matter EoS . Up to nuclear saturation den- 
sity, the nuclear matter EoS can be experimentally con- 
strained. If one wants to extrapolate it to densities above 



II. CRITERION FOR STABLE HYBRID STARS 

A. Connected hybrid branch 

A compact star will be stable as long as the mass M of 
the star is an increasing function of the central pressure 
Pcent [E]. There will therefore be a stable hybrid star 
branch in the M(R) relation, connected to the neutron 
star branch, if the mass of the star continues to increase 
with p C ent when the quark matter core first appears, at 

Vcent — Pcrit j 



connected hybrid dM 
star branch dp ccnt 



> 



pi 



(2) 



This derivative is discontinuous at p C rit; we evaluate it for 
central pressures infinitesimally above the critical pres- 
sure. In this region, the contribution of the quark mat- 
ter core is dominated by the energy density discontinuity 
Ae, since the quark core is not large enough for the slope 
( c qm) of e(p) to have much influence on the mass and 
radius of the hybrid star. 
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FIG. 2: Mass-radius relations for hybrid stars made of HLPS nuclear matter [2j with the transition to quark matter at n cr it = 4no 
(no = 0.16 fm -3 is nuclear saturation density) and various values of the energy density discontinuity Ae. Thick dashed lines are 
for stable hybrid stars, thin dashed lines show unstable branches. Case (a): Cq M = 1/3; the stable hybrid branch disappears as 
Ae rises beyond Ae s tab- Case (b): Cq M = 1; a disconnected stable hybrid branch is present for Ae > Ae s tab- 



If Ae is small then quark matter has a similar energy 
density to that of nuclear matter, and we expect hybrid 
stars to be stable when the quark matter core first ap- 
pears. If Ae is too large, then the star becomes unstable 
as soon as the quark matter core appears, because the 
pressure of the quark matter is unable to counteract the 
additional downward force from the gravitational attrac- 
tion that the additional energy in the core exerts on the 
rest of the star. We therefore expect that a stable hybrid 
star branch will exist if Ae is less than a threshold value 
Ae s tat>; 



QM EoS that gives 
connected hybrid branch 



Ae < Ae 



stab • 



(3) 



the mass will start rising again when the central pressure 
rises sufficiently far above p CT it , producing a disconnected 
branch of stable hybrid stars. This constitutes a "third 
family" |20l |2"T| of compact stars besides neutron stars 
and white dwarfs. Such a feature has been noticed in 
M(R) calculations for specific quark matter models, for 
example kaon condensed stars [52], quark matter cores 
from perturbative QCD [23 , and color superconducting 
quark matter cores (23] • Also Ref. (35] demonstrates the 
transition from connected branch to disconnected branch 
via twin star collapse. Here we will study the generic 
features of a quark matter EoS that give rise to this phe- 
nomenon. 



The situation is illustrated in terms of mass-radius plots 
in Fig. [2^ a), which contains M(R) curves for the HLPS 
nuclear matter EoS [5] with a transition to quark matter 
at ?icrit = 4no . In this case, we assume quark matter with 
Cq M = 1/3 and we vary Ae. The hybrid branch becomes 
unstable (mass decreasing with central pressure) when 
Ae is greater than a critical value Ae s t a b which in this 
case is 0.727e cr it. 



B. Disconnected hybrid branch 



In Fig. |2jb), we illustrate the possibility of a second, 
disconnected, branch of stable hybrid stars at Ae > 
Ae s tab- The only difference from Fig.[2ja) is that we have 
increased the speed of sound in quark matter to Cq M = 1. 
The critical value of Ae for a connected hybrid branch is 
of course unchanged, but if Ae is above that value then 



C. Calculating the critical energy density 
discontinuity Ae sta b 

Applying the criterion in Eq. ^ for arbitrary transi- 
tion pressures, we obtain a function Ae sta b(« C rit) which 
depends only on the nuclear matter equation of state. 
To evaluate it for a given nuclear matter EoS, we choose 
a value of n cr it and scan through values of the energy 
density discontinuity Ae to find the largest value that 
gives a stable hybrid star when the central pressure is 
pushed just above p cr it- We then vary n cr ; t to obtain 

Ae sta b(«crit)- 

In Fig. [3] we show the results of a scan through val- 
ues of Ae for n cr it = 1.5 uq. In this calculation we fix 
Cq M = 1/3, but recall that the value of this parameter 
does not affect Ae s t a b- In the example displayed in Fig. [3] 
the mass M of the star increases with central energy den- 
sity for small Ae, but as Ae increases the rate of increase 
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FIG. 3: Calculation of Ae sta b: Mass versus the central energy 
density for quark matter equations of state with various energy 
density discontinuities Ae. The nuclear matter EoS is NL3 and 
the transition density is fixed at 1.5 no. The hybrid star branch 
becomes unstable when Ae is greater than Ae ata b ~ 0.6e C rit, 
because the mass then decreases with central pressure when the 
quark matter core first appears. 



of M with ^central becomes smaller, until Ae reaches the 
value Ae s tab ~ 0.6e cr it at which point the mass becomes 
a decreasing function of central energy density when the 
quark matter core appears, indicating an unstable star. 
There is then no stable hybrid star branch connected to 
the nuclear star M(R) curve. However, as seen in the 
lowest line in Fig. [3] (for Ae = 0.7e cr ; t ), even if the mass 
at first decreases with central pressure, it may then turn 
around and start rising again, indicating the presence of a 
stable hybrid branch disconnected from the nuclear mat- 
ter branch. This is the scenario illustrated in Fig. [2|b) . 



D. "Phase diagram" for hybrid stars 

In Fig. [4] we plot a "phase diagram" for hybrid stars, 
where the control parameters are the baryon density n cr j t 
of nuclear matter at the transition to quark matter and 
the discontinuity Ae in the energy density at the transi- 
tion. This diagram is obtained by repeating the calcula- 
tion of Fig. [3] for many values of n cl .; t . We present results 
for a softer nuclear matter EoS (HLPS) and a stiffer one 
(NL3). 

The region below the solid curve (Ae < Ae s tab( ri crit)) 
is where there is a stable hybrid star branch connected to 
the hadronic star branch. We see that, apart from critical 
densities close to the upper limit n max where the hadronic 
star itself becomes unstable, Ae sta b is always about 60% 
to 80% of the energy density e cr j t of the nuclear matter at 
the hadronic-quark transition. (For HLPS we also mark 
the density n — 5.458 no at which HLPS becomes un- 
physical because its speed of sound goes above 1.) Nat- 
urally, as n cr it comes close to the value n max at which 



the hadronic star itself is unstable, the window for a con- 
nected branch of stable hybrid stars becomes small and 
at n cr it = n m ax it vanishes. 

Above the solid line in Fig. [3] is the region Ae > 
Aegtab^crit) where there is no connected hybrid star 
branch. The almost-vertical dashed lines in that region 
demarcate, for two different values of Cq M , the regions 
where there is a separate hybrid star branch, discon- 
nected from the hadronic branch (see Fig.[2jD). The prop- 
erties of such stars depend on the behavior of the quark 
matter EoS above the deconfining transition, so we show 
results for two values of Cq M . 

We see in Fig. [4] that the disconnected branch occurs 
more readily when the transition density is low, and 
the speed of sound in quark matter is high. It is also 
slightly favored by a lower value of Ae. Such features 
were noticed in the context of stars with mixed phases 
in Ref. [TU], which pointed out that the disconnected 
branch can be understood as follows. When the quark 
matter core first appears, its greater density creates a 
strong gravitational pull on the nuclear mantle, which 
the pressure of the core cannot support, so the star is 
unstable. However, if the pressure of the core rises fast 
enough with energy density, a larger core with a slightly 
higher central density may be able to sustain the weight 
of the mantle above it. This explains why disconnected 
branches occur more readily when the transition density 
is low: the nuclear mantle is then less dense, and easier 
for the quark core to support. It also explains why the 
vertical line separating the "disconnected hybrid branch" 
region from the "no hybrid branch" region moves to the 
right as Cq M increases, since c 2 = dp/de is a measure of 
stiffness: if Cq M is larger then the pressure of the core 
grows more quickly as its density grows, making a large 
core better able to support the nuclear mantle. Finally, 
it explains why the line has a slight negative slope: larger 
Ae makes the quark core heavier, increasing its pull on 
the nuclear mantle, and making the hybrid star more un- 
stable against collapse. 

Ref. |10j . assuming a mixed phase, conjectured that 
the third family (disconnected branch) exists when the 
speed of sound in quark matter is less than that in the 
mixed phase. In the case we study, with no mixed phase, 
the relevant quantity would be the difference between 
the speed of sound in quark matter Cq M and the speed 
of sound c NM in nuclear matter at the phase transition. 
However, by comparing Fig. [4] with Fig. [5] we see that, 
although c NM rises with ro cr ; t , it is not the case that 
the disconnected hybrid branch disappears exactly when 
c nm > c qm- F° r example, for NL3 nuclear matter the 
speed of sound never reaches 1, but when NL3 nuclear 
matter is combined with quark matter with Cq M = 1 the 
disconnected hybrid branch still disappears at a finite 
transition density. 
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FIG. 4: "Phase diagram" for hybrid star branches in the mass-radius relation. The control parameters are the nuclear matter 
density n cr i t at the transition to quark matter, and the energy density discontinuity Ae at the transition. When Ae < Ae sta b 
there is a connected hybrid branch. This is independent of the the speed of sound Cq M in the quark matter. When Ae > Ae sta b 
there may be a disconnected hybrid branch if the critical density is below a threshold value which depends on Cq M . We show 
results for a softer nuclear matter EoS (HLPS j2|) and a harder one (NL3 |16j). 
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FIG. 5: Speed of sound squared c 2 in the HLPS and NL3 nu- 
clear matter equations of state. At n < 4no HLPS is a softer 
EoS with a lower speed of sound. At higher density it becomes 
harder and then acausal. 



III. MAXIMUM MASS OF HYBRID STARS 

A. Maximum mass and central energy density 

The generic quark matter EoS (fij) involves three pa- 
rameters, p cr j t (or equivalently n cr ;t7, Ae, and Cq M . For 
each value of these parameters that gives stable hybrid 
stars, we calculated the maximum mass M max of the hy- 



brid star. The results of Ref. |26[ [27] lead us to expect 
that, for a given nuclear matter EoS, the maximum mass 
is mostly determined by the central energy density of the 
heaviest star and the speed of sound in the central regions 
of that star. Specifically, 

M max = y(c cont )£^/ t 2 . (4) 

The function y(c cen t) can be obtained from Ref. |H] 
(their Eq. (24) and the associated table). To test Eq. Q, 
we follow Ref. [26 , 27 and plot the maximum mass M max 
that a stable star can have, as a function of the central 
energy density in that star e ccn t- We use the NL3 nu- 
clear matter EoS and repeat the calculation for a range 
of quark matter equations of state, varying the critical 
density n cr it for the transition to quark matter, the quark 
matter speed of sound Cq M , and the energy density dis- 
continuity Ae at the transition. The range of allowed 
values of Ae, giving stable hybrid stars, is from zero to 

Ae sta b(«crit)- 

The results are shown in Fig. [7J There are solid curves 
showing M max (e ccnt ) according to Eq. Q for c 2 ccnt = 1, 
0.6, and 1/3. Nuclear matter equations of state at high 
density have c 2 close to 1, hence the maximum masses for 
pure nuclear matter stars (we show the result for APR 
[T] nuclear matter as well as NL3 and HLPS) lie close to 

the C cc„t = 1 line - 

The triangles, squares, and circles show the masses 
of hybrid stars for a quark matter EoS with relatively 
low transition densities: n crit = 1.5 n (open symbols) 
or n cr it = 2 no (solid symbols) , and for each of these 
we vary Cq M and Ae. The low transition density means 
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(a) HLPS + QM (b) NL3 + QM 



FIG. 6: Contour plot of the mass of the heaviest hybrid star as a function of quark matter EoS parameters n cr i t , Ae, Cq M , for 
HLPS (left panel) and NL3 (right panel) nuclear matter. The thin (red), medium (green) and thick (blue) lines are for a nuclear 
to quark transition at n cr it = 1.5no, 2no, and 4no, respectively. The dashed lines are for stars on a disconnected branch. 
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FIG. 7: Mass of the heaviest hybrid star as a function of its 
central energy density, for various quark matter equations of 
state Q. The curves are predictions of Ref. for stars whose 
core-region speed of sound squared is 1, 0.6, and 1/3. Pure 
nuclear matter stars for the NL3 and APR equations of state 
are also plotted. 



that these hybrid stars have large quark matter cores, 
so the the speed of sound in the central part of the star 
is Cq M , so the maximum masses should fall on the lines 
given by Eq. (ffl) with c cont = c QM . We see that on the 
whole this is the case. However, for 7i cr ; t = 2.0 n Q and 
Cq M = 1/3 (solid triangles) the points fall slightly be- 
low the predicted Cq M = 1/3 line. This implies that in 
these stars the effective speed of sound in the core is even 



lower than l/\/3. However, the speed of sound in the nu- 
clear matter mantle at n ~ 2no is greater than (see 
Fig. [5]). It seems reasonable to argue that the first-order 
phase transition, at which dp/de — c 2 vanishes, acts like 
a "soft" region where c is small [12], and drags down the 
average value of c con t . All hybrid stars have such a tran- 
sition region, but in these stars the quark matter cores 
are smaller than in the other cases, so the nuclear-quark 
matter boundary is closer to the center of the star, and 
we conjecture that it plays a more important role in de- 
termining the maximum mass. 



B. How heavy can a hybrid star be? 

Using the very generic pameterization (JT|) of the quark 
matter EoS, it is possible to get hybrid stars that are 
heavy enough to be consistent with the recent measure- 
ment of a star of mass 2 Mq (28]. In Fig. [6] we show con- 
tour plots for the the maximum masses of hybrid stars 
as a function of the parameters of the quark matter EoS, 
for both the HLPS and NL3 nuclear equations of state. 
The solid lines are contours of given maximum mass for 
hybrid stars on a connected branch. These naturally end 
at Ae = Ae s tab (from Fig. [4J. The dashed lines are for 
stars on a disconnected branch. 

In both cases we see that, as expected from Fig. [4] 
hybrid stars will be heavier if the energy density discon- 
tinuity Ae is smaller (so the gravitational pull of the core 
does not destabilize the star) and the speed of sound in 
quark matter is higher (so the core is stiffer and can sup- 
port a heavy star) . 

For the softer HLPS EoS, whose pure nuclear matter 
star has a maximum mass of 2.15 Af (Table hybrid 
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stars can be heavier than the pure nuclear star. This 
occurs when the transition density is low, and the quark 
core is a large fraction of the star. In the M(R) relation, 
the hybrid branch can have a positive dM / dR, so its 
M(R) curve looks similar to that of a pure quark matter 
star, rising to a maximum mass star which is heavier and 
larger than the heaviest pure HLPS star. For the harder 
NL3 EoS, whose pure nuclear matter star has a maximum 
mass of 2.77 M Q , hybrid stars are always lighter then the 
heaviest nuclear matter star. 



C. The quark matter mass fraction 

For hybrid stars, it is natural to ask is how much of 
the mass of the star consists of quark matter. We de- 
fine the quark matter mass fraction f q = M core /Mst ar . 
In Fig. |6ja), along the thick (blue, n cr i t = 4rio) contours 
of constant M star , fq 

varies from about 60% at low Cq M 
and Ae to about 70% at high Cq M and Ae. If the transi- 
tion pressure is lower then we expect the quark fraction 
to be larger: along the medium-thickness (green) mass 
contours for n cr ; t = 2n , f q varies from about 90% at low 
Cq M and Ae to about 96% at high Cq M and Ae. 

In Fig.[6jh), along the medium-thickness (green) mass 
coutours for n cr ; t = 2n Ql f q varies from about 50% at 
low Cq M and Ae to about 80% at high Cq M and Ae. 
Again, if the transition pressure is lower then the quark 
fraction is larger: along the thin (red) mass contours for 
"■crit = l-5ri0) fq varies from about 80% at low Cq M and 
Ae to over 90% at high Cq M and Ae. 



IV. CONCLUSIONS 

We studied hybrid stars where there is a sharp inter- 
face between two phases with different equations of state. 
We called the two phases "nuclear matter" and "quark 
matter", but our conclusions are valid for any first-order 
phase transition between two phases with different en- 
ergy densities. 

We found that in this situation there is a simple cri- 
terion for determining whether there exists a stable con- 
nected branch of hybrid stars. For a given nuclear matter 
EoS, the only relevant features of the quark matter EoS 
are the transition density/pressure and the discontinuity 
in energy density Ae at the transition. A stable con- 
nected branch of hybrid stars exists if Ae is less than a 
critical value Ae s t a b which is a function of the nuclear 
matter EoS and the density at which the transition oc- 
curs, and can be easily calculated for any nuclear EoS. 
Thus for each nuclear matter EoS there is a "phase di- 
agram" (e.g. Fig. |4| in which the line Ae s t a b(^crit) de- 
lineates the region where there is a connected hybrid 
branch. In the examples we studied Ae s t a b is typically 
around 60% to 80% of the nuclear matter energy density 
at the transition, with a rapid drop to zero when the crit- 



ical density becomes close to the maximal central density 
of a pure nuclear matter star. 

To study properties of hybrid stars beyond the central 
pressure where they first appear, we used a fairly generic 
classical-ideal-gas-type parameterization of the quark 
matter EoS at densities beyond the transition, and came 
to the following conclusions. 

(1) Even if there is no connected hybrid star branch, 
there may be a disconnected one if the transition 
density is low enough so that the nuclear mantle can 
be supported by a large enough quark matter core (see 
Fig. [J. 

(2) We confirmed that the relationship Q between the 
maximum mass of a star and its central energy density, 
which was proposed in Ref. [26 , holds for hybrid stars 
(see Fig. [7). 

(3) We found that it is reasonably easy to get heavy 
hybrid stars (more than 2M 9 ) for reasonable parame- 
ters of the quark matter EoS. It requires a not-too-high 
transition density (n cr ;t ~ 2no), low enough energy den- 
sity discontinuity Ae < 0.5e cr it, and high enough speed 
of sound Cq M > 0.4). A value of Cq M that is well above 
1/3 is an indication that the quark matter is strongly 
coupled. For a hard nuclear matter EoS such as NL3 it 
was somewhat easier to make heavy hybrid stars. For 
details see Fig. [6] This confirms what was seen in previ- 
ous calculations for specific models of quark matter (e.g., 
[ 1 1 1 IT5HT5] ) : in the absence of theoretical or experimental 
constraints on the quark matter EoS, one can fairly eas- 
ily vary the unknown parameters of that EoS to obtain 
heavy hybrid stars. For now, it seems clear that without 
theoretical advances that constrain the form of the quark 
matter EoS, measurements of gross features of the M(R) 
curve such as the maximum mass will not be able to rule 
out the presence of quark matter in neutron stars. 

The generic parameterization of the quark matter EoS 
that we use in this work provides a general framework for 
comparison and empirical testing of models of the quark 
matter EoS. Any particular model can be characterized, 
at least at densities close to the transition, in terms of 
the parameters ncrit j Ae, and Cq M , and its predictions 
for connected hybrid stars will follow from its position 
in the phase diagram (Fig. |4j). If the form of the nu- 
clear matter EoS were established then measurements of 
the M(R) relation of neutron stars could be directly ex- 
pressed as constraints on the values of our quark matter 
EoS parameters. 

Our discussion of stars with central pressure well above 
the transition density relied on the assumption of a 
density-independent speed of sound, which is a useful 
starting point for general comparisons of quark matter 
models, as well as providing specific examples of quark 
matter equations of state that can yield heavy hybrid 
stars (Fig. [6|. If observations of M(R) for heavy stars 
turned out to be inconsistent with our parameteriza- 
tion, or if theorists were able to show that the speed 
of sound in quark matter has significant density depen- 
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dence, then our approach could be further generalized 
(at the penalty of introducing more parameters) to allow 
for that, and to allow for mixed or percolated phases as 
noted in Refs. [TUl fTT]. 
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Note that we have introduced an additional parameter 
A with mass dimension 3-/3. The value of A can 
be varied without affecting the energy-pressure relation 
(All. When constructing a first-order transition from 



some low-pressure EoS to a high-pressure EoS of the form 
(All, we must choose A so that the pressure is a mono- 



tonically increasing function of fig (i.e. so that the jump 
in tib at the transition is not negative). The derivative 
with respect to fis yields 



tib(hb) = (l + /3)A(i B 



(A3) 



and using p ~ hbtt-b — £, we obtain the energy density 



Appendix A: Ideal gas equation of state 

Here we briefly recapitulate (see, e.g., Ref. [T^]) the 
construction of a thermodynamically consistent equation 
of state of the form in Eq. (JT|) 



e (p) = £o + \p ■ 



(Al) 



e(ji B ) =B + pA f i 1 B 



(A4) 



Then Eq. ( A2 1 gives energy density as a function of pres- 
sure 



We start by writing the pressure in terms of the chemical 
potential 



Mb (p) 



P + B 
A 



(A2) 



e(p) = (1 + 0)B + /3p 



(A5) 



which is equivalent to Eq. ( Al ) with 1/c 2 = /? and £o 
(1 + 0)B. 
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